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Abstract. If gamma ray bursts are highly coUimated, radiating into only a small 
fraction of the sky, the energy requirements of each event may be reduced by several 
(up to 4-6) orders of magnitude, and the event rate increased correspondingly. The 
large Lorentz factors (F > 100) inferred from GRB spectra imply relativistic beaming 
of the gamma rays into an angle ^ 1/T. We are at present ignorant of whether there 
are ejecta outside this narrow cone. 

Afterglows allow empirical tests of whether GRBs are well-coUimated jets or spher- 
ical fireballs. The bulk Lorentz factor decreases and radiation is beamed into an ever 
increasing solid angle as the burst remnant expands. It follows that if gamma ray 
bursts are highly coUimated, many more optical and radio transients should be ob- 
served without associated gamma rays than with them. In addition, a burst whose 
ejecta are beamed into angle Cm undergoes a qualitative change in evolution when 
rCm ^ 1: Before this, F oc r~^^^, while afterwards, F c>c exp{—r/r^). This change 
results in a potentially observable break in the afterglow light curve. 

Successful application of either test would eliminate the largest remaining uncer- 
tainty in the energy requirements and space density of gamma ray bursters. 



The ejecta from gamma ray bursts must be highly relativistic to explain the 
spectral properties of the emergent radiation [1,4]. The gamma rays we observe 
are therefore only those from material moving within angle l/F of the line of sight, 
and offer no straightforward way of determining whether the bursts are isotropic 
emitters or are beamed into a small angle. (Here F is the bulk Lorentz factor of 
expansion.) 

Afterglow emission at longer wavelengths is expected to arise later in the evo- 
lution of the burst than the original gamma rays. It therefore offers at least two 
ways of testing the burst beaming hypothesis. 
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Burst and Afterglow Event Rates 



First, because T is lower at the time of afterglow emission than during the GRB 
itself, the afterglow cannot be as coUimated as the GRB can. This implies that the 
afterglow event rate should exceed the GRB event rate substantially if bursts are 
strongly beamed. Allowing for finite detection thresholds, 
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where A^i, N2 are the measured event rates above our detection thresholds at our 
two frequencies; N12 is the rate of events above threshold at both frequencies; and 
fli, Q2 are the solid angles into which emission is beamed at the two frequencies. 

A full derivation of this result and discussion of its application is given in [6]. 
Rather than reproduce it, I will refer the reader to that paper and will here discuss 
the second test more fully than was possible in [6] . 

Dynamical Calculations: Numerical Integrations 

The second test is based on differences between the dynamical evolution of 
beamed and isotropic bursts. We explore the effects of beaming on burst evo- 
lution using the notation of [5]. Let Fq and Mq be the initial Lorentz factor and 

ejecta mass, and Cm the opening angle into which the cjccta move. The burst energy 
is E = FoMoC^Cm/^j where we assume a unipolar jet geometry. Let r be the radial 
coordinate in the burster frame; t, tco, and the time from the event measured in 
the burster frame, comoving ejecta frame, and terrestrial observer's frame; and / 
the ratio of swept up mass to Mq. 

The key assumptions in our beamed burst model arc that (1) the energy and 
mass per unit solid angle are constant at angles 9 < Cm from the jet axis and 
zero for 9 > Cm (see [2] for an alternative model); (2) the energy in the ejecta is 
approximately conserved; (3) the ambient medium has uniform density; and (4) the 
cloud of ejecta + swcpt-up material expands in its comoving frame at the sound 
speed Cs = cj \/3 appropriate for rclativistic matter. The last of these assumptions 
implies that the working surface of the expanding remnant has a transverse size 
~ Cm'^ + Cs^co- The evolution of the burst changes when the second term dominates 
over the first. 

The full equations describing the burst remnant's evolution are then 
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/ = — / r'^lm{r)p{r)dT , (2) 
Mq Jo 

= T^{Cm + Cstco/Ctf ^ n{Cm + tco/V3tf , (3) 



T = {To + f) /Jl + 2Tof + f ^ jTo/2f , (4) 



t = r/c , tco= tdt'/T , and = t dt' /2V'^ . (5) 
Jo Jo 

These equations can be solved by numerical integration to yield /(r), r(r), and 
t^{r). Figure 1 shows r(r) from such integrations for an illustrative pair of models 
(one beamed, one isotropic). 
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FIGURE 1. Dependence of the bulk Lorentz factor T on the burst expansion radius for an 

isotropic burst and a burst beamed into an opening angle Cm = 0.01 radian. Both bursts follow 
a r oc r""^/^ evolution initially, but the beamed burst changes its behavior at F « 100 « l/Cmi 
beyond which its Lorentz factor decays exponentially with radius. 

The emergent synchrotron radiation can also be calculated if we assume an elec- 
tron energy spectrum and assume that electrons and magnetic fields have constant 
fractions of the equipartition energy density. For illustrative purposes, we again 
follow the assumptions in [5]. The electron energy spectrum is N{£) oc S"^, i.e. a 
power law with equal energy per decade, so that the synchrotron spectrum peaks 
where r = 0.35, rising as z/^^^ at low (optically thick) frequencies and falling as 
at high (optically thin) frequencies [3]. The relevant equations are a straight- 
forward modification of equations 11-20 of [5]. Figure 2 shows the peak fiux density 
as a function of observed frequency for the models used in figure 1. We caution the 
reader that more recent electron energy spectra grounded in observations (e.g. [7]) 
may be more reliable. We hope to incorporate such spectra in our calculations in 
future. 

Dynamical Calculations: Analytic Integrations 

The most interesting dynamical change introduced by beaming is a transition 
from a power law F oc r~^/^ to an exponentially decaying regime F oc exp(— r/r^,). 
This can be derived by considering the approximate evolution equations for the 
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FIGURE 2. The dependence of the peak fliix density fi, on observed frequency u for the same 

pair of bursts. The electron spectrum follows the model of Paczynski & Rhoads (1993). The 
peak in the synchrotron emission for this model occurs at the frequency where optical depth 
effects become important. The predicted break in the power law caused by beaming should be 
observable. Similar breaks occur in the dependence of and Vpeak with time, and are expected 
to be a generic feature of beamed GRB afterglow models. 



regime where (a) I/Tq < f < Fq, so that F ^Jf^/QJ] and (b) Cgtco > CmT 
(corresponding to / > 9roC^): 



df/dr^^cit'ip , dtJdr^J^ , dtjdr^^ . (6) 
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It follows that 



This is easily integrated to obtain 



The constant of integration becomes negligible once Cgtco ^ Cm^i so that equa- 
tion 8 becomes / oc t^^. It is then clear from equations 6 that /, F, t^o, and 
will all behave exponentially with r in this regime. Retaining the constants of 
proportionality, we find 



/ oc exp(2r/rp) where 
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Further algebra yields T oc exp(— r/rp) and oc / oc exp(2r/rp), so that F oc ^0 . 
Thus, while the evolution of F(r) changes from a power law to an exponential at 
F ~ 1/Cm, the evolution of t©(r) changes similarly. The net result is that r{t^) has 
a power law form in both regimes, but with a break in the slope from F oc f^^^^ 
when F > to F oc t^^^"^ when F < 1/Cm- 

Of course, F is not directly observable, and we ultimately want to predict observ- 
ables like the frequency of peak emission Um, the flux density F^^^n at u^, and the 
angular size 9 of the afterglow. With the electron energy spectrum described above, 
the relevant power law scalings before beaming becomes dynamically important are 

T^m ~ f^^^^, F^^rn ~ ^©^^^^, and 9 ~ t%^. At late times, Urn ~ t©S -Pl,m ~ t®^^^, and 
1/2 

9 ^ . Our numerical integrations confirm these relations, though the transition 
between the two regimes is quite gradual for u^. 

Combining these scalings with the spectral shape yields predictions for the light 
curve at fixed observed frequency. The most dramatic feature is in the light curve 
shape for u > Um, which changes from Fj,^© ~ © ~ t^^. These exponents 

are generally sensitive to the assumed electron energy distribution in the blast wave. 



Conclusions 

Establishing whether or not gamma ray bursts are beamed will be valuable in 

understanding source populations and burst mechanisms. There are at least two 
potentially observable consequences of beaming. 

(1) The event rate for afterglows should exceed that for bursts substantially if 
bursts are strongly beamed. A quantitative comparison of rates at two frequencies 
yields quantitative limits on the ratio of beaming angles. 

(2) The dynamical evolution of a beamed burst remnant changes qualitatively 
when F < 1/Cm- The resulting changes in the light curves could be observed. 
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